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Abstract 

In this paper, analytical expressions correlating the volumetric flow rate to 
the pressure drop are derived for the flow of Carreau and Cross fluids through 
straight rigid circular uniform pipes and long thin slits. The derivation is based 
on the application of Weissenberg-Rabinowitsch-Mooney-Schofield method to 
obtain flow solutions for generalized Newtonian fluids through pipes and our 
adaptation of this method to the flow through slits. The derived expressions 
are validated by comparing their solutions to the solutions obtained from 
direct numerical integration. They are also validated by comparison to the 
solutions obtained from the variational method which we proposed previously. 
In all the investigated cases, the three methods agree very well. The agreement 
with the variational method also lends more support to this method and to 
the variational principle which the method is based upon. 

Keywords: fluid mechanics; rheology; non-Newtonian fluids; Carreau; Cross; 
pipe; slit; Weissenberg-Rabinowitsch-Mooney-Schofield equation; variational 
method. 


1 Introduction 

There are many fluid models that have been developed and employed in the recent 
decades to describe and predict the bulk and in situ rheology of non-Newtonian 
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fluids. Amongst these, Carreau, and to a certain degree Cross, are distinguished by 
their popularity and widespread use especially in modeling the rheological behavior 
of biological fluids and polymeric liquids. Both are four-parameter models that 
depend on the low-shear and high-shear viscosities, a characteristic time or strain 
rate constant and a flow behavior exponent. They are usually used to describe the 
time-independent shear-thinning category of the non-Newtonian fluids. 

These widely used models provide good match to the experimental data in many 
flow situations, such as the flow in the blood arteries and through porous media as 
well as rheometric measurements, and hence they are popular in various biological, 
technological and industrial disciplines, such as biosciences and engineering, reservoir 
engineering and food processing. For example they are systematically used, in their 
different variants and forms and their modified versions, to model and predict the 
flow of biological fluids like blood, and polymeric liquids like Xanthan gum and 
polyacrylamide gel solutions [1-22], 

To the best of our knowledge, no analytical solutions that correlate the volumetric 
flow rate to the pressure drop in confined geometries, specifically pipes and slits, 
have been reported in the literature for these fluids [6, 23-25] despite their wide 
application. This is mainly due to the rather complicated expressions that result 
from applying the traditional methods of fluid mechanics to these models in any 
analytical derivation treatment. Therefore, the users of these models either employ 
empirical approximations or use numerical approaches which normally utilize mesh- 
based techniques like finite element and finite difference methods. 

In this paper, we make an attempt to derive fully analytical solutions for the 
flow of these fluids through straight rigid circular uniform tubes and thin long slits. 
We use a method attributed to Weissenberg, Rabinowitsch, Mooney and Schofield 
[26], and may be others, and hence we call it WRMS method. The method, as 
reported in the literature, is customized to the flow in circular pipes; and hence we 


2 


adapt it to the flow in thin slits to obtain flow relations for this type of conduit 
geometry as well. 

In fact there are two objectives to the present paper. The first is the derivation 
of the analytical expressions, as outlined in the previous paragraph, which is useful, 
and may even be necessary, in various rheological and fluid mechanical applications. 
The second, which is not less important, is the solidification and support to our 
recently proposed [27-29] variational method which is based on optimizing the total 
stress in the flow conduit by applying the Euler-Lagrange principle to find totally 
or partly analytical flow relations for generalized Newtonian fluids through various 
conduit geometries. As we will see, the results of the newly derived formulae in the 
present paper for the flow of Carreau and Cross fluids through pipes and slits agree 
very well with the results obtained from the variational method. Since the two 
methods are totally independent and are based on completely different theoretical 
and mathematical infrastructures, they provide support and validation to each 
other. 

The plan for this paper is that in section § 2 we present the WRMS method for 
the flow of generalized Newtonian fluids through pipes and derive its adaptation for 
the flow through slits. The WRMS is then applied in section § 3 to derive analytical 
relations for the flow of Carreau and Cross fluids through pipes, while its adaptation 
is used in section § 4 to derive these relations for the flow through slits. The derived 
analytical expressions are then validated in section § 5 by numerical integration 
and by comparison to the flow solutions which are obtained from the variational 
method. The paper is ended in section § 6 with short discussion and conclusions 
about the purpose and the achieved objectives of this study and possible future 
extensions to other models. 
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2 Method 


First, we should state our assumptions about the flow, fluid and geometry and 
mechanical properties of the employed conduits. In this investigation we assume a 
laminar, incompressible, isothermal, steady, pressure-driven, fully-developed flow of 
a time-in dependent, purely-viscous fluid that is properly described by the generalized 
Newtonian fluid model where the viscosity depends only on the contemporary rate 
of strain and hence it has no deformation-dependent memory. As indicated earlier, 
the pipe is assumed to be straight rigid with a uniform and circularly shaped cross 
sectional area while the slit is assumed to be straight rigid long and thin with a 
uniform cross section. It is also presumed that the entry and exit edge effects and 
external body forces, like gravity, are negligible. As for the boundary conditions, we 
assume no-slip at the tube and slit walls [30] with the flow velocity profile having a 
stationary derivative point at the symmetry center line of the tube and symmetry 
center plane of the slit which means that the profile has a blunt rounded vertex. 

We first present the derivation of the general formula for the volumetric flow 
rate of generalized Newtonian fluids through pipes that satisfy the above-stated 
assumptions, as depicted in Figure 1. The derivation is attributed to Weissenberg, 
Rabinowitsch, Mooney, and Schofield and hence we label it with WRMS. The 
outline of this derivation is a modified version of what is in [26] which we reproduce 
here for the purpose of availability, clarity and completeness. We then follow this by 
adapting the WRMS method from the pipe geometry to the long thin slit geometry 
which we also call it WRMS method. The difference between the two will be obvious 
from the context. 

The differential volumetric flow rate in a differential annulus between r and 
r + dr, as depicted in Figure 1, is given by 

dQ = v2n rdr (1) 
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Figure 1: A diagram depicting the geometry of the straight rigid circular uniform 
pipe that is used in the current study. 

where Q is the volumetric flow rate, r is the radius and v = v(r) is the fluid velocity 
at r in the axial direction. The total flow rate is then obtained from integrating the 
differential flow rate between r = 0 and r = R where R is the tube radius, i.e. 



On integrating by parts we obtain 



( 3 ) 


The first term on the right hand side of the last equation is zero because at the lower 
limit r = 0 , and at the upper limit v = 0 due to the no-slip boundary condition 
[30]; moreover 


dv = — 7 r dr 


( 4 ) 
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where y r = 7 (r) is the rate of shear strain at r. On eliminating the first term, 
substituting for dv and changing the limits of integration, Equation 3 becomes 


f R 

Q = tt r 2 ^ r dr 

Jo 


( 5 ) 


Now, by definition, the shear stress at the tube wall, t w , is the ratio of the force 
normal to the tube cross section, F±, to the area of the internal surface parallel to 
this force, An, that is 


. . F± 7rR 2 Ap RAp 

r^r( r = fl)^r R =- = 1 - 5r = — 


( 6 ) 


where L is the tube length, and A p is the pressure drop across the tube. Similarly 


r(r) = r r = 


rAp 

~2L 


(0 < r < R) 


( 7 ) 


Hence 


T r r 

Tr R 


2 7? 2 T r 2 

r = — 


T 


( 8 ) 


R 


and 


dr R R 

-— = — dr = — dr r 

dr r Tr Tr 


( 9 ) 


On substituting r 2 and dr from the last two equations into Equation 5 and changing 
the limits of integration we get 


Q = IT 


rTR R 2 t 2 R , nR 3 

2 d T r — 


CTr 


'o r R T R 


T 


T 2 -f r dT r 


( 10 ) 


R JO 


that is 
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( 11 ) 


Q 


7 r 


R 3 



where it is understood that 7 = 7 ,. = 7 (r) and r = r r = r(r). 

Next, we derive a general formula for the volumetric flow rate of generalized 
Newtonian fluids through a long thin slit, depicted in Figure 2, by adapting the 
WRMS method, which is described in the last paragraph, as follow. The differential 
flow rate through a differential strip along the slit width is given by 


dQ = vWdz 


( 12 ) 


where W is the slit width and v = v(z) is the fluid velocity at 0 in the x direction 
according to the coordinates system demonstrated in Figure 2. Hence 


Q. 

w 


f +B 


vdz 


(13) 


l-B 


where B is the slit half height with the slit being positioned symmetrically with 
respect to the plane z = 0. On integrating by parts we get 


Q r 1 +B 

W = [vz] ~ B ‘ 


nv + B 


zdv 


(14) 


'V-B 


The first term on the right hand side of the last equation is zero due to the no-slip 
boundary condition at the slit walls, and hence we have 


Q 

w 


t* v +B 


zdv 


(15) 


'V_s 


Now, if we follow a similar argument to that of pipe, the shear stress at the slit 
walls, t w , will be given by 


bn — hi-B 


An 


2BW Ap BAp 
2WL = L 


(16) 
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Figure 2: A diagram depicting the geometry of the straight rigid uniform long thin 
slit that is used in the current study. 


where L is the slit length and A p is the pressure drop across the slit. Similarly we 
have 


r, = 


zAp 


(17) 


where 7% is the shear stress at z. Hence 


z Bt z 

T z = ~f;T± B ->• 2 = - 

B t ±b 


—y dz — 


Bdr z 

t±b 


(18) 


We also have 


dv 


Bdr z 


lz = -~r dv = -j z dz = -7*- 

dz t±b 


(19) 


Now due to the symmetry with respect to the plane z = 0 we have 


t b = t +b = r_ B (20) 

On substituting from Equations 18 and 19 into Equation 15, considering the 
symmetry with respect to the center plane z = 0 and changing the limits of 
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integration we obtain 


that is 


Q_ 

W 


rr + s B 


't-b 


t±b 


~lz 


Bdr z 

r±B 



( 21 ) 



where it is understood that 7 = = 7 (z) and r = r z = r{z). 


( 22 ) 


3 Pipe Flow 

In the following two subsections we apply the WRMS method to derive analytical 
expressions for the flow of Carreau and Cross fluids in straight rigid circular uniform 
pipes. 

3.1 Carreau 

For Carreau fluids, the viscosity is given by [31, 32] 

fi = -=m + {fi 0 - /h) (1 + aV) (n_1)/2 (23) 

where /i 0 is the low-shear viscosity, \i % is the high-shear viscosity, A is a characteristic 
time constant, and n is the flow behavior index. This, for the sake of compactness, 
can be written as 


M = - = ft + ,5(l + AV)” 72 (24) 
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where 5 = (/x 0 — /ij), and n' = (n — 1). Therefore 


9 





r = 7 


(25) 


ft+ 5(1 +A V )” 72 

From WRMS method we have 


u=r^ dr p6) 

If we label the integral on the right hand side of Equation 26 with / and substitute 
r from Equation 25 into / we obtain 



Now, from Equation 25 we have 


A,+5 (i+Avr 72 


dr 


(27) 


dr = 


ft + <5(1 + AY )” 72 + n'ttV (1 + AV) 


(n' — 2)/2 




(28) 


On substituting from Equation 28 into Equation 27 and changing the integration 
limits we get 


f7iu 


I = 


7 


2 „, 2 \"V 2 


di + $ (l + -^y 2 ) 


2 r 


ft + a (1 + A 2 7 2 )“ 72 + n'5AV (1 + A 2 7 2 ) ( ”'- 2)/2 

(29) 


G?7 


where 7 ^ is the rate of shear strain at the tube wall. On solving this integral 
equation analytically and evaluating it at its two limits we obtain 
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i 2 _,2 \3n'/2 


/ = 


<r [3A 4 (3 n' 2 + 5n' + 2) y 4 - 3n'A 2 7 2 + 2] (1 + A^y 2 ) 
3A 4 (9 n' 2 + 18 n' + 8 ) 

I^P [A 4 ( 2 n' 2 + 5ra' + 3) p - n'\ 2 p w + 1] (1 + A 2 y 2 )” 


2A 4 (n' + 1) in' + 2 ) 


lij5 [A 4 (n /2 + 5 n' + 6 ) y 4 - n'\ 2 p w + 2 ] (1 + A 2 y 2 )" + 


A 4 ( n' + 2 ) (n' + 4) 

25 3 iip 


+ 


+ 


4 


3A 4 (9n /2 + 18 n' + 8 ) 2 A 4 (n' + 1 ) (n' + 2 ) A 4 (V + 2 ) (n' + 4) 


30 ) 


For any given set of fluid parameters, the only unknown that is needed to compute 
/ from Equation 30 is r y w . Now, by definition, through the application of the main 
rheological equation to the flow at the tube wall, we have 


7~w 


(31) 


that is 


hi + <5 (1 + ^ 2 lw) 


n '/2 




RAp 

2 L 


(32) 


From the last equation, 7 ^ can be obtained numerically by a simple numerical 
solver, based for example on a bisection method, and hence / is computed from 
Equation 30. The volumetric flow rate is then obtained from 


Q = 


nR 3 I 


(33) 


which is fully analytical solution apart from computing the value of j w . However, 
since obtaining 7 ^ numerically can be easily achieved to any required level of 
accuracy, as it only depends on very simple and reliable solution schemes like the 
bisection methods, this does not affect the analytical nature of the solution and its 
accuracy. 
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3.2 Cross 


For Cross fluids, the viscosity is given by [33] 


(1 

7 


/^o 


— Ah + 


(34) 


where /i Q is the low-shear viscosity, /i l is the high-shear viscosity, A is a characteristic 
time constant, m is an indicia! parameter, and 6 = (/i G — Therefore 


t = 7 /h + 


(35) 


1 + X m Y' 

If we follow a similar procedure to that of Carreau by applying the WRMS method 
and labeling the right hand side integral with / we get 


rTw f A 

„3 ' 0 


1 J 0 7 l /il + 1 + \ m r 


dr 


Now, from Equation 35 we have 


(36) 


dr = ^ + 


m5 \ m Y 


d'f 


(37) 


1 T A m y m (i _)_ A m 7 m ) 2 ; 

On substituting from Equation 37 into Equation 36 and changing the integration 
limits we get 


/ 




1 + X m Y‘ 


+ 1 + A m 7 m 


m<5A m 7 m A 

(1 + A m 7 m f ) 


c?7 


(38) 


On solving this integral analytically and evaluating it at its two limits we obtain 
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/ = 


{2<r [—m (2 f 2 + 5/ + 3) + 4 g 2 + 2m 2 ] + 12 m5 2 /iig (m — r/) + 12 m 2 6g 2 g 2 + 3m 2 g 3 g 3 } 7 , 


v 4 

IW 


12m 2 g 3 

{<5 3 (m 2 - 6 m + 8 ) + 3m<J 2 /^ (m - 4) + 3 m 2 5g 2 } 2 F 1 (l, +; -/) 7 


12 m 2 


(39) 


where 


/ = A m 7 ™ , 9 = 1 +f (40) 

and 2 -Pi is the hypergeometric function of the given argument with its real part 
being used in this evaluation. As before, we have 


that is 




(41) 



* 1 

1 + A m 7™ y 




RAp 

2L 


(42) 


From the last equation, 7 ^ can be obtained numerically, e.g. by a bisection method, 
and hence / is evaluated. The volumetric flow rate is then computed from 


Q = 


7 tR 3 I 


(43) 


which is fully analytical solution, as explained in the Carreau case. 


4 Slit Flow 

In the following two subsections we apply the adapted WRMS method to derive 
analytical relations for the flow of Carreau and Cross fluids in straight rigid uniform 
long thin slits. 
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4.1 Carreau 


For Carreau fluids, the viscosity and shear stress are given by Equations 24 and 25 
respectively. Now, from the adapted WRMS method for slits we have 


Qt, 


'yrdr 


(44) 


2 WB 2 

If we label the integral on the right hand side of Equation 44 with /, substitute for 
r and dr from Equations 25 and 28 into /, and change the integration limits we 
obtain 




1=1 7 Z U + *(1 + AY)” 72 Li + 6 (1 + a 2 7 2 )" 72 + (1 + aV) ( "'- 2)/2 

Jo L - 1 L 

(45) 

On solving this integral analytically and evaluating it at its two limits we obtain 


d'y 


I = 


+ 

+ 


n'5 2 ^ w [ 2 F x Q, 1 - n'\ |; -A 2 7 2 ) - 2 F\ (A, -n'; §; -A 2 7 2 )] 

A 2 

(1 + n') S 2 'jl 2 (§ j -n'; -A 2 7 2 ) 


n'Sfa[ 2 F, (|, 1 - f; !; -A 2 7 2 ) - 2 *i (±,-£; |; -A 2 7 2 )] 

A 2 

(2 + n') 5/i i7 3 2 -Pi (|, -y; |;-A 2 7 2 ) + 


(46) 


where 2 -E 1 is the hypergeometric function of the given arguments with the real part 
being used in this evaluation. Now, from applying the rheological equation at the 
slit wall we have 


fii + 5 (1 + A 2 7 2 ) 


n '/2 


w 


B A p 


(47) 
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From the last equation, can be obtained numerically by a simple numerical solver 
and hence / is computed. The volumetric flow rate is then obtained from 


Q 


2 WB 2 I 


(48) 


4.2 Cross 

If we follow a similar procedure to that of Carreau flow in slits by applying the 
adapted WRMS method and labeling the right hand side integral with / we get 


/ = 


7 [Hi + 


dr 


(49) 


1 + \ m y‘ 

where dr is given by Equation 37. On substituting from Equation 37 into Equation 
49 and changing the integration limits we get 


r r w 2 ( s 
~ J 0 7 v' + 1 + \ m r 


Hi + 


m6 X m Y 


1 + A m 7 m (l + 


d'y (50) 


On solving this integral equation analytically and evaluating it at its two limits we 
obtain 


I = 


[3h 2 (■ m-g ) - {5 2 (m - 3) + 2 mh/ 7 } g 2 2 C (l, 1 + ^; -/) + 6 mS^g + 2 mtfg 2 ] 7 , 


6 mg 2 


(51) 


where 


/ = A m 7^ , g = 1 + / (52) 

and 2 F\ is the hypergeometric function of the given argument with its real part 
being used in this evaluation. As before, from applying the rheological equation at 
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the wall we have 



* ) 
1 + A m 7™ J 


'Yw 


B Ap 
L 


(53) 


From this equation, 7 ^ can be obtained numerically and hence / is computed. 
Finally, the volumetric flow rate is obtained from 


Q 


2 WB 2 I 


(54) 


5 Validation 

We validate the derived equations in the last two sections for the flow of Carreau and 
Cross fluids through pipes and slits by two means. First, we validate the derived / 
expressions (i.e. Equations 30, 39, 46 and 51), which are the main source of potential 
errors in these derivations, by numerical integration of their corresponding integrals 
(i.e. Equations 29, 38, 45 and 50 respectively). On comparing the analytical solutions 
to those obtained from the numerical integration of /, we obtained virtually identical 
results in all the investigated cases. This validation eliminates the possibility of a 
formal error in the derivation of the analytical expressions of / but does not provide 
a proper validation to the basic WRMS method since the numerical integration does 
not eliminate possible errors in the fundamental assumptions and basic principles 
upon which the WRMS method is based and the derivation steps that lead to the 
integrals / and subsequently to Q. The second way of validation, which will be 
explained in the following paragraphs, should provide this sort of validation as it 
is based on comparing the WRMS solutions to solutions obtained from a totally 
different method, namely the variational method, which was already validated. 

In [27] a variational method based on applying the Euler-Lagrange variational 
principle to find analytical and semi-analytical solutions for the flow of generalized 
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Newtonian fluids in straight rigid circular uniform pipes was proposed. There, the 
method was applied to the Carreau and Cross fluids, among other types of fluid, 
where a mixed analytical-numerical method was developed and used to find the flow 
of these two fluids in pipes. Later on [28, 29], the variational method was further 
validated by extending it to the slit geometry and to more types of non-Newtonian 
fluids, namely Ree-Eyring and Casson. Elaborate details about all these issues, as 
well as other related issues, can be found in the above references. 

Extensive comparisons between the variational and the WRMS analytical so¬ 
lutions have been carried out as part of the current investigation to validate the 
analytical solutions on one hand and to add more support to the variational method 
on the other. In all the investigated cases, which vary over a wide range of fluid 
and conduit parameters, very good agreement was obtained between the variational 
and the WRMS analytical solutions. In Figures 3-6 we present a representative 
sample of the investigated cases where we compare the WRMS analytical solutions 
with the variational solutions for the flow of Carreau and Cross fluids through 
pipe and slit geometries. As seen in these figures, the agreement between the two 
methods is excellent which is typical of all the investigated cases. The main source 
of error and departure between the two methods is the heavy use of numerical 
bisection solvers and numerical integration in the implementation of the variational 
method. The presence of the hypergeometric function in both the WRMS and 
variational solutions is another source of error since this complicated function may 
not converge satisfactorily in some cases when evaluated numerically, and hence 
becomes problematic and a major source of error. It is worth mentioning that the 
solutions obtained from the numerical integration of I integrals are not shown on 
Figures 3-6 because they are virtually identical to the analytical solutions. 
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Figure 3: Comparing the WRMS analytical solutions to the variational solutions for 
the flow of Carreau fluids in straight rigid circular uniform tubes with (a) n = 0.65, 
/i G = 0.1 Pa.s, Hi = 0.005 Pa.s, A = 1.5 s, L = 0.85 nr and R = 0.09 nr; and (b) 
n = 0.9, Ho = 0.08 Pa.s, Hi — 0.001 Pa.s, A = 2.0 s, L — 0.5 nr and R = 0.02 nr. In 
both sub-figures, the vertical axis represents the volumetric flow rate, Q, in m 3 .s -1 
while the horizontal axis represents the pressure drop, A p, in Pa. The average 
percentage relative difference between the analytical and variational solutions for 
these cases are about 0.37% and 0.60% respectively. 
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Figure 4: Comparing the WRMS analytical solutions to the variational solutions for 
the flow of Cross fluids in straight rigid circular uniform tubes with (a) m = 0.83, 
Ho = 0.22 Pa.s, Hi — 0.033 Pa.s, A = 6.65 s, L = 0.95 m and R = 0.008 m; and (b) 
m = 0.5, Ho = 0.15 Pa.s, Hi — 0.009 Pa.s, A = 7.9 s, L = 1.65 m and R = 0.027 m. 
In both sub-figures, the vertical axis represents the volumetric flow rate, Q, in m 3 .s _1 
while the horizontal axis represents the pressure drop, A p, in Pa. The average 
percentage relative difference between the analytical and variational solutions for 
these cases are about 0.55% and 0.32% respectively. 
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(a) (b) 


Figure 5: Comparing the WRMS analytical solutions to the variational solutions 
for the flow of Carreau fluids in straight rigid uniform long thin slits with (a) 
n = 0.75, Ho = 0.17 Pa.s, p t = 0.009 Pa.s, A = 2.5 s, L = 1.3 m, W = 1.0 nr and 
B = 0.012 m; and (b) n = 1.0, p a = 0.09 Pa.s, pi = 0.012 Pa.s, A = 3.0 s, L = 1.1 m, 
W = 1.0 m and B = 0.015 m. In both sub-figures, the vertical axis represents the 
volumetric flow rate, Q, in m 3 .s -1 while the horizontal axis represents the pressure 
drop, Ap, in Pa. The average percentage relative difference between the analytical 
and variational solutions for these cases are about 0.70% and 0.06% respectively. 




Figure 6: Comparing the WRMS analytical solutions to the variational solutions for 
the flow of Cross fluids in straight rigid uniform long thin slits with (a) m = 0.45, 
p Q = 0.08 Pa.s, pi = 0.003 Pa.s, A = 0.75 s, L = 0.75 nr, W = 1.0 nr and 
B = 0.005 nr; and (b) m = 0.75, p Q = 0.03 Pa.s, pi = 0.005 Pa.s, A = 5.0 s, 
L = 1.25 nr, W = 1.0 nr and B = 0.015 nr. In both sub-figures, the vertical axis 
represents the volumetric flow rate, Q, in m 3 .s _1 while the horizontal axis represents 
the pressure drop, A p, in Pa. The average percentage relative difference between 
the analytical and variational solutions for these cases are about 0.89% and 0.40% 
respectively. 
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6 Conclusions 


In this paper, a method based on the application of the Weissenberg-Rabinowitsch- 
Mooney-Schofield (WRMS) equation for the flow of generalized Newtonian fluids in 
straight rigid circular uniform pipes and its extension to straight rigid long thin 
uniform slits has been developed and used to find analytical flow solutions for 
Carreau and Cross fluids which do not have known analytical solutions in those 
geometries. The main analytical expressions were verified by numerical integration 
to rule out potential errors in the formality of derivation. The analytical solutions 
were then thoroughly compared to the solutions of another method, namely the 
variational approach which is based on the application of the Euler-Lagrange 
principle and hence it is totally different and independent from the WRMS method. 
Excellent agreement was obtained in all the investigated cases. 

The first thing that has been achieved in this investigation is obtaining fully 
analytical solutions for the flow of Carreau and Cross fluids in the above mentioned 
geometries. These solutions provide a better alternative to the use of numerical 
techniques, which are currently the only available means, since the analytical 
solutions are easier to obtain, less prone to error and highly accurate. The second 
is that these analytical solutions provide more support to the previously proposed 
variational method for obtaining solutions for the flow of generalized Newtonian 
fluids in confined geometries. In fact this agreement should serve as mutual 
validation since these methods are totally independent and are based on very 
different theoretical and mathematical infrastructures and hence they lend support 
to each other. 

It should be remarked that the derivation of the pipe and slit analytical relations 
for the Carreau model can be easily extended to the more general Carreau-Yasuda 
model by replacing ‘2’ with ‘a’ in the exponents of the Carreau constitutive relation. 
The WRMS method may also be applied to other non-Newtonian fluid models which 
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do not have analytical solutions. However, in some cases the definite integrals, /, 
may require numerical evaluation by using one of the numerical integration methods 
such as quadratures, Simpson, trapezium or midpoint integration rules if fully 
analytical solutions are very difficult or impossible to obtain. Such semi-analytical 
solutions will not only be much easier to obtain and less prone to errors than their 
numerical counterparts, but can also be as accurate as any potential analytical 
solutions for all practical purposes due to the availability of very reliable numerical 
integration schemes that can be employed. These schemes are easier to implement 
and more reliable than the alternative numerical techniques like those based on 
discretization schemes. 
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Nomenclature 


7 shear rate (s _1 ) 

shear rate at conduit wall (s^ 1 ) 

8 Ho~ Hi (Pa.s) 

A characteristic time constant (s) 

H dynamic shear viscosity of generalized Newtonian fluid (Pa.s) 

Hi high-shear viscosity (Pa.s) 

Ho low-shear viscosity (Pa.s) 

r shear stress (Pa) 

tb shear stress at slit wall (Pa) 

Tr shear stress at tube wall (Pa) 

t w shear stress at conduit wall (Pa) 

A\\ area of conduit internal surface (m 2 ) 

B slit half height (m) 

/ 

F± force normal to conduit cross section (N) 

2 -F] hypergeometric function 

9 1 + / 

/ definite integral expression (Pa 3 .s -1 for pipe and Pa 2 .s _1 for slit) 
L length of conduit (m) 

m indicia! parameter in Cross model 

n flow behavior index in Carreau model 

n' n — 1 

A p pressure drop (Pa) 
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Q volumetric flow rate (m 3 .s -1 ) 
r radius (nr) 

R tube radius (m) 

v fluid velocity in the flow direction (m.s -1 ) 
W slit width (m) 

z coordinate of slit smallest dimension (m) 
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